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We present a self-consistent method for deriving moment equations for kinetic models of polymer 
dynamics. The Doi model |lj of liquid-crystalline polymers with the Onsager excluded-volume 
potential is considered as an example. To lowest order, this method amounts to a simple effective 
potential different from the Maier-Saupe form. Analytical results are presented which indicate that 
this effective potential provides a better approximation to the Onsager potential than the Maier- 
Saupe potential. Corrections to the effective potential are obtained. 
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I. INTRODUCTION 

Kinetic theory is a powerful analytical tool for describing the dynamics of dilute and semi-dilute solutions of 
polymers in terms of a diffusion equation for the particle distribution function or, equivalently, by the full system of 
moment equations. In general, the moment system has to be truncated at some level. The problem of the "closure 
approximation" is very well-known in the literature, and an enormous amount of suggestions have been analyzed in 
the case where each moment couples only to a few higher-order moments. However, for some kinetic equations, the 
time evolution of each moment couples to an infinite set of higher moments, so that further analytical work is often 
precluded since closure approximations are less studied for this case. In this paper we present a simple method that 
allows to derive moment equations with a finite coupling valid for a wide class of kinetic equations. 

In order to be specific, we consider a particularly important example: the Doi theory of liquid-crystalline polymers 
(LCP), subject to the Onsager excluded-volume potential [Q. As it is well-known, in this model each moment 
equation depends on an infinite set of higher-order moments. In the original work this problem was treated in 
two steps: First, the Onsager potential was replaced by a different, phenomenological potential of the Maier-Saupe 
type which gives raise to a coupling to the next higher moment only. In the second step, the "decoupling" 
approximation was used to solve the resulting closure problem for the second moment. Subsequent extensive studies 
were focused on improvements of the second step 1^-^. At the same time, we are not aware of improvements on 
the first step and closure approximations are limited to the Maier-Saupe potential up to now. However, it would be 
desirable to deal with the true Onsager potential, not only because it becomes exact in the limit of low concentrations 
of perfectly rigid rod-like molecules, but also because it contains no phenomenological parameters and therefore gives 
more quantitative predictions. In addition, the Onsager potential is preferred in the study of the influence of flow on 
the isotropic-nematic transition, since it gives a clear-cut prediction of the range of coexistence of the equilibrium 
isotropic and nematic phase, both in stationary |Q| and non-stationary jsj flows. The method, which we propose 
in this work, leads to an approximation of the Onsager potential, which, to the lowest order, is at the same time 
as simple as the Maier-Saupe potential but also closer to the true Onsager potential. Moreover, corrections to this 
approximation can be obtained in a systematic manner. 

II. THE DOI MODEL 

Let ^{u; t) be the probability distribution function for a rigid rod-like polymer molecule to be oriented parallel to 
the unit vector u. The time evolution of in the presence of flow and the Onsager excluded-volume potential was 
given by Doi 0| and may be written as: 



dti! ^-n-[ux{K.- uii,)] + n ■ D.ipn ij^ij^ ) • (i) 
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Here TZ = ux d/du is the rotational operator, d/du the gradient on the unit sphere, k the gradient of the velocity, D,- 
the rotational diffusivity, 5 / the functional derivative and A = -\- Ai is the free energy functional per molecule 
divided by k^T, 

Ao=\nv-l + {\nip{u)) (2a) 
U 



A, = ^{Wl-{u-wY)). (2b) 

U — 2bL'^v is the reduced excluded- volume, 2h and L are the diameter and the length of the rod-like polymeric 
molecules, respectively, and v is the number of molecules per unit volume. Here and below we use the following 
notations for averages: (/(«)) = / f{u)ijj{u)du, and {{f{u,w))) —JJ f{u,w)tp{u)ip{w)dudw. Aq describes the loss 
of entropy with molecular alignment, while Ai expresses the Onsager free energy of steric interaction in the second 
virial approximation B . Following Doi and Edwards pG] , the rotational diffusivity is approximated by 



(3) 



where D^, the rotational diffusion coefficient for a rod in an isotropic, semi-dilute solution of like rods, is related 
to the rotational diffusion constant for a dilute solution, -DrOj by Z?r = cDyo{i'L)~^ with an empirical coefficient c. 
Nonlinearity of Eq. (|l|) in ip brought about by the potential (^) reflects the mean-field nature of the Onsager theory 
of the excluded-volume effect. The self-consistent potential, identified by Doi is related to the free energy of 
interaction, V{u) — kBT6Ai/6ip{u). Various phases of the LCP are conveniently described by the order parameter 
S = {uu — (1/3)1), where 1 is the unit tensor. It is reasonable therefore to look for approximate formulations of the 
dynamics in terms of the order parameter alone. However, as mentioned above, the time evolution equation for S 
couples to an infinite number of moments of tp- In the derivation given by Doi, this difficulty was circumvented by 
replacing the Onsager potential (|2b| ) by a different, phenomenological expression of the Maier-Saupe type j^: 

AfS = „^ _ : S, (4) 

where gq and oi are parameters independent of "0. A further separate treatment of the diffusivity (H) is also necessary. 
A compact presentation of the entire development is given by Doi and Edwards [0 . The Doi model with the Maier- 
Saupe potential (Q) constitutes the basic kinetic model of LCP used by many authors for analytical studies to derive 
equations for the order parameter. As is well-known, the kinetic equation (|^) with the potential (Q) does not give a 
closed equation for the order parameter but contains also the higher-order moment (uuuu) , and therefore constitutes 
a further problem of closure. The original Doi approach was based on the decoupling approximation for the fourth- 
order moments of ip in terms of S. Improvements on the decoupling approximation are currently under active research 



III. GENERATING MOMENT EQUATIONS 

In this communication, we demonstrate that a different self-consistent treatment of the kinetic equation (|l|) is 
possible. Modifications concern only the relaxational part of the Eq. (|l|), specifically, the excluded-volume potential 
( pb| ) and the diffusivity (||), and therefore we consider the case k = in the sequel to simplify notations. Specifically, 
we employ the cumulant expansion of the potential ( pb[ ) and the diffusivity (p|). The leading term of this expansion 
results in an effective potential that differs from the Maier-Saupe potential (H)^ and which contains a non-polynomial 
dependence on the order parameter S. 

In the second virial approximation, the free energy of interaction, Ai can be written as Ai — {v /2){{(3{u,w))) . If 
only excluded-volume interactions are present, the second virial coefficient (3 corresponding to the Onsager expression 
( pb| ) is of the form f]{u,w) = (3{{u-w)'^), with /3{x) = 2hl?\J\ — x. Specifically, expanding /3(a;) in a Taylor series and 
interchanging summation and averaging in this expansion, we get ((/3(x))) = X]n=o ^"((■^"))' '^l^^re a„ are numerical 
coefficients. Each average ((a;")) can be represented in terms of cumulants {{x^Jjc of order k < n. Resummation of the 
series leads to 

((/3(x))) =/?(((:.))) + ^-p^J /3(2™)(((^))) + ..., (5) 
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where is the 2m-th. derivative of (3 and ellipses denote terms including third or higher-order cumulants as 

factors. Therefore, the functional Ai (|2b| ) can be split as Ai = A^^"^ + B, where A'"^^ is the free energy, corresponding 
to the total neglect of second and higher-order cumulants in each term of the expansion, 

4'^ = |v/l (6) 

In terms of the order parameter S, A\ may be rewritten as A^^' = {U/V6)y/1 - (3/2)3 : S. By the mean value 
theorem, it is easy to see that A'"^^ gives an upper bound to Ai, Ai < A\^\ for the present case of excluded-volume 
interactions. The functional B contains the higher-order cumulants. While all powers of the second cumulants are 
displayed in Eq. (|^), in general it is not a priori clear whether it is more important to keep powers of the second 
cumulant or higher cumulants. However, we generally expect the linear term in the second order cumulant to be 
most important. The corresponding term m = 1 in Eq. (|^), — {l/2)((x'^))cl3" {((x))), gives the first non-vanishing 
contribution to B, 

4" = ■ - ■■ (1 - ■■ {uu))-"^ . (7) 

Keeping only the first N cumulants in the expansion (^, the functional (^) is approximated by non-linear functions 
of the first 2N moments of ijj. Inserting this approximation in the time evolution equation (|^) amounts to an 
approximation of the self-consistent Onsagcr potential V in terms of a polynomial of order 2N in u but with non- 
linear coefficients. In this approximation, the time evolution of the 2n-th moment contains only the first (27i + 2iV)-th 
moments. With this, moment equations can be generated, that approach the original equations in a systematic way, 
thereby containing only a finite number of moments at each stage. 



IV. TESTING THE APPROXIMATION 

Clearly, the above procedure is most valuable if the first terms, A\^^ , etc. , already provide a good approximation to 
the full expression Ai. While the general validity of ^4^^^ as a good approximation to Ai is a rather delicate problem, 
it should be mentioned that it is so at least in two limiting cases. Namely, for the isotropic state, the value of A']^'' 
differs from Ai for less than 5%, while in the fully ordered state the approximation (^) becomes exact. Moreover, on 
the submanifold of distribution functions of the form 

a 

ipa{u) = — — cosh(aM • n), (8) 

Att smh a 

where n is an arbitrary unit vector, and < a < cxd, the functional ^^^^ turns out to approximate Ai very well for all 
values of the parameter a between the isotropic state, a = 0, and the fully ordered state, a = oo. To show this, we 
plot in Fig. 1 the functions Ai{a), A\^\a) and A^^^a) + A^^\a), that result upon inserting the ansatz (^) into (p^), 
®; (0); respectively. Note, that Ai{a), A^j^\a) and ^^'(q!) can be calculated analytically. For convenience, we plot 
the functions against the scalar order parameter, defined as S* = ^/ (3/2)3 : S. Including A^^-* does not only reduce the 
error of the approximate value of Ai in the isotropic state to 1.5%, but improves the accuracy of the approximation 
over the whole range of 5*. For comparison, we included in Fig. 1 also the free energy A^^, corresponding to the 
Maier-Saupe expression (^, thereby choosing the undetermined constant so that the limit of the fully ordered state 
is matched correctly. Note, however, that in any case A^^ decays asymptotically like 1/a, for a ^ 1, whereas Ai 
and A'"^^ behave like l/^/a in this regime. We included in Fig. 1 also the derivative of the above functions, since 
they are related to the self-consistent potential V. Fig. 1 shows that also the derivative of A^^^ provides a good 
approximation to the derivative of Ai, with correct limiting behavior near the isotropic and fully ordered state. Note 

(2) 

that including the first correction A\ yields excellent agreement to the true Onsager prediction. The Maier-Saupe 
potential captures the main features but, besides an undetermined constant, shows the wrong behavior near the fully 
ordered state. The ansatz (H) , originally proposed by Onsager , is known to approximate the equilibrium distribution 

very well. Therefore we conclude that A^^"^ represents a good approximation to Ai, at least on a representative subset 
of distribution functions. 
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V. THERMODYNAMIC CONSISTENCY 



It is worth mentioning that the presentation given so far can easily be cast into the recently developed GENERIC 
formalism of nonequilibrium thermodynamics ||ll|,|l^. In the absence of potential forces, the example of rigid dumb- 
bells, which are equivalent to the model of rigid rods, is formulated within the GENERIC formalism in Ref. The 
mean field potentials considered above can be included in a straightforward manner, if we recognize that Ai — —Si, 
where Si is the entropic contribution per molecule to the free energy of interaction divided by fee- Formulating the 
original model as well as the approximations within the GENERIC formalism guarantees that our treatment is in 
accordance with the principles of nonequilibrium thermodynamics. This becomes especially important if the present 
model is considered in nonisothermal situations. For example, the structure of the GENERIC formalism requires the 
polymeric contribution to the elastic stress to be 

S S 

(^al3 = "ivkBTSali ~ vk^T {{u X 'R — )aUp). (9) 

Eq. (0) agrees with the result of Doi pOj, obtained upon varying the free energy functional. 



VI. THE LOWEST ORDER APPROXIMATION 

In the sequel, we will adopt the lowest order approximation A = A^^^ = + ^i'^'' > where is given by Eq, 
(^, and Aq is given by Eq. (|2^). This amounts to neglect of all higher order correlations in Eq. (|^), or, equivalently 
setting i? = 0. Substituting instead of A into Eq. (|l]), we derive 

Uuu : {uu) 



TZtp - tpTZ 



2-\/l — (uu) : (uu) 
It is now possible to identify the self-consistent potential as 



(10) 



FW(n)=(^)-i=^^=i^, (11) 

which can be compared to the expression obtained from inserting the Maier-Saupe free energy (|^) into Eq. (|l|) 

Vms (u) = a2 — aiUkBTuu : (uu), (12) 

where 02 is an arbitrary constant. The normalized equilibrium solutions to the Eq. ([l0| ) are ipcli = 
Z-iexp[-F(iVfcBT]. 

The rotational diffusivity (^) is related to the free energy of interaction, since D, = D, [^Ai/{U/2)] . Substituting 
Ai = a''^^ gives 

d[^'' = (37rV32)A-[l - (3/2)5 : S]-\ (13) 
The diffusion coefficient I?!^'' (E3|) is positive in the entire physically meaningful range of the order parameter S. 



Expression (|l^ ) should be compared with the Doi phcnomenological result: 

Ad- A[l- (3/2)5: 5] (14) 

While we have not found an argument which of the two powers, —1 or —2, is more consistent, it should be stressed 
that our derivation of the diffusion coefficient does not need any further assumptions or adjustable parameters, while 
the derivation of Eq. (|l|) §0 requires the matching of D,-, resp. A^^^ in both, the isotropic and the fully ordered 

state. Due to its relation to a\^\ the diffusion coefffcient Z?|^'' (^3|) has a correct limit in the fully ordered state 

{Di./d[^^ = as soon as 5 : 5 = 2/3 in the ordered state), while the opposite limit of the isotropic state {D^ — A) 
is matched within 8%. Again, the first correction (^ reduces the error in this limit to less than 3%. 

If we adopt (^ and approximate by D^^^ (p^, the time evolution of the order parameter S can be derived from 
Eq. (O) by the so-called Prager procedure 
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dtS = -Gd'-^^ S + 6d[^'' , ^ =(S-(uu) - S : (uuuu)), (15) 
' Vl- (3/2)5: 5^ \ / \ \ ) 

with U' = U/VQ- This expression differs from the result of Doi 0| not only in the diffusion coefficient and in the 
reduced excluded-volume U due to the undetermined constant in the Maier-Saupe potential (p^, but contains a 
non-polynomial dependence on the order parameter 5, which becomes important in the nematic state. 



VII. CONCLUSION 



We have presented a systematic procedure that allows to derive approximate moment equations for the Doi model of 
LCP, which contain only a finite number of higher order moments. The first approximation for the Onsager excluded- 
volume interaction results in an effective potential ( p+l ) proportional to uu, but different from the Maier-Saupe form 
( p^ ) and without free parameters. We find indications, that ( pi] ) approximates the true Onsager potential better than 
the Maier-Saupe potential. For higher accuracy, the first correction seems to be the most important contribution. 
All these approximations are in accordance with nonequilibrium thermodynamics. 

Note, that we have not addressed the problem of solving the resulting kinetic equations or "closing" the moment 



equations. This work is currently under preparation. Nevertheless, for comparing Eq. (15) to the corresponding 
equation with the Maier-Saupe potential, we follow Refs. |l|Jlo[] and consider the decoupling approximation S : 
(uuuu) — S : {uu){uu). If the order parameter is assumed to be of the form Saf3 — S{t)[nanp — Sap], the relaxation 
equation for the scalar parameter 5' is found to be as follows: 

dtS = -6D^. — , (16a) 



where U' = U/VG. Due to a non-polynomial character of the function A^-^^^ (161), the relaxation equation (16a) differs 
formally from the Landau-de Gennes counterpart derived by Doi for the Maier-Saupe potential. Expansion of the 
function ( [l6b| ) around 5 = reproduces the result of Doi for the Maier-Saupe potential, subject to a rcnormalization 
of the strength of the excluded-volume potential, and a difference in the coefficient in front of the 5^ term. Moreover, 
the relaxation implied by A^^^ ( |l6bD is qualitatively similar to the one given by Doi result and distinguishes the same 
three regimes. For U <Ui, A^^has only one minimum at 5 = 0, so that the system finally becomes isotropic. For 
Ui < U < U2, Si second local minimum occurs. The system either becomes isotropic or nematic depending on the 
initial value of S. Finally, for U > U2, the isotropic state becomes unstable and the system always approaches a 
nematic state. 

Due to the undetermined constant oi in the Maier-Saupe potential ([T^), the Doi theory predicts the values Ui 
and U2 also in terms of ai. On the contrary, the self-consistent potential ( pT| ) contains no free parameters, so that 
Ui = S^/'^VS/iVS - 1) « 6.22 and U2 = 3^6 « 7.34 may directly be compared to the values in the Onsager theory 
Ui = 8.38 and U2 = 10.67. However, the values of the order parameter Si — 1/4 and 5*2 — 1/2 at Ui resp. U2 in the 
Doi theory do not depend on ai and can therefore be compared to the values predicted by (pT|): Si = (\/3— 1)/2 « 0.37 
and 52 = l/\/2 « 0.71. 

In Fig. 2, the equilibrium order parameter 5cq is shown as a function of iy/i'2- The lower solid line shows the 
prediction of the Doi theory, whereas the upper solid line corresponds to the approximation ([ll|). For the Maier- 
Saupe potential, 5cq 7^ is given by 



For the free energy ( |l6q ), Scq is given implicitly as the solution to the algebraic equation 

1 + 5cq - 25^q ^ 1^2 ^^^^ 



For large values of v/v2, the solution ( [Tq ) approaches the value 5oq = 1 and asymptotically behaves like the solution 
of the Doi theory 5eq ~ l—V2/v, for large u. Note, that the decoupling approximation corrects the asymptotic behavior 
of the Maier-Saupe potential near the fully ordered state. 
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As is well-known, the detailed form of the interaction potential can have significant effect on the behavior of the 
order parameter in the nematic phase . Specifically, the amount of order at the transition is known to be much 
smaller in the Maier-Saupe theory than in the Onsager model. For comparison, we included in Fig. 2 the values of the 
order parameter obtained from minimizing the true Onsager free energy numerically [l5| ], where vi now corresponds 
to the true nematic transition. Although the analysis of the phase transitions via the dynamical approach is affected 
by the use of the decoupling approximation, the prediction of the self-consistent approach is much closer to the true 
Onsager values than is the Maier-Saupe potential. 

Finally, it should be mentioned that approximations to the Onsager potential like Eqs. ^ and (Q) can also be used 
in the case of potential flows, following the approach of Thirumalai without additional assumptions. 

To summarize, we have developed a direct approach to the Doi model with the Onsager potential. We have 
demonstrated that the resulting kinetic equation has much in common with the Doi model with the phenomenological 
Maier-Saupe potential. Corrections to the approximation developed here can be found in a systematic way from Eqs. 

to (0) by taking into account higher order correlations. The approach to derive self-consistent moment equations 
is applicable to other kinetic equations which can be cast into the form (|l|). 
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FIGURE CAPTIONS 



Fig. 1: Free energy of excluded-volume interaction for distribution functions (g) plotted against the scalar order 
parameter ^(q;). From top to bottom: approximation v4^^'' (|6|), with first correction A''^^ + A^^' true 
Onsager expression Ai (2b) and the Maier-Saupe free energy A^^ (^), when the limits 5 = and 5 = 1 are 
matched. In the inset, the derivative of the above functions is shown as a function of 5(a). The order of the 
curves from top to bottom is the same. 

Fig. 2: The equilibrium order parameter Scq as a function of i'/i'2- The figure shows the behavior due to the Maier- 
Saupe potential (p7[), lower curve, and the solution of ( p^ ) corresponding to A^-^\ upper curve, in the decoupling 
approximation. Full circles indicate the order parameter for the true Onsager potential in the static case (from 

PI)- 
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